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SUMMARY 


A rapid method for estimating the downwash Behind swept-ving air- 
planes is presented. The Basic assumption is that of a flat horizontal 
sheet of vortices trailing Behind the wing. The integrations for the 
downwash are handled in a manner similar to Both Multhopp*s and 
Weissinger's approximate integrations in their span-loading calculations. 
The principal effects of rolling -up of the wake are treated as correc- 
tions to the flat-sheet wake. A simple approximate correction for the 
effect of the fuselage is applied. The agreement with availaBle experi- 
mental data taken Behind airplane models is good. Computing forms are 
included together with charts of pertinent functions, so as to enable 
simple direct application. 


INTRODUCTION 


The downwash induced By a lifting wing has, in the past. Been pre- 
dicted By considering the wing as a lifting line with a vortex sheet 
trailing aft of the wing in a horizontal plane. It was assumed that 
spanwise distribution of vorticity did not change with downstream posi- 
tion and that the sheet did not roll up Behind the wing. With these 
assumptions, a procedure for detennining downwash is given in refer- 
ences 1 and 2. In references 1 and 2, the wing span loading is approxi- 
mated By several horseshoe vortices. The total downwash is the sum of 
the downwashes of the horseshoe vortices. It is apparent that such a 
procedure can Be extended to swept wings By using swept horseshoe vor- 
tices. The arithmetic of this procedure is, however, rather tedious 
and laborious. In reference 3^ a more rapid method in the form of an 
inf luence -coefficient approach is presented for the downwash at the 
center of the wake. References 1 and 2 also investigated the limitations 
of representing the lifting surface By a lifting line, and of the effects 
of the rolling -up of the trailing sheet. It was concluded that Both 
effects were negligible for the then conventional airplane configurations. 

At the present time, the use of low-aspect-ratio plan forms and 
occasionally of further rearward positions of the tail has made neces- 
sary a re-examination of the assumption that the trailing vortex sheet 
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can be considered nonrolling-up. An analysis of the rolling-up process 
is given in reference 4 which reveals that the trailing sheet becomes 
rolled-up at shorter distances behind the wing as (l) aspect ratio 
decreases, (2) lift coefficient increases, and (3) span loading increases 
outboard and decreases inboard. It is apparent that the downwash fields 
determined on the assumption of the flat trailing. vortex sheet or a com- 
pletely rolled-up sheet (the simplified cases) omit wings of aspect 
ratio of about two to four at moderate or high ' s . 

The purposes of this report are, (l) to make available an influence - 
coefficient type of method of computing the downwash behind swept wings 
having arbitrary spanwise loading, a procedure that will be quicker and 
simpler to use than methods summing up the downwash due to elemental 
horseshoe vortices, (2) to estimate the principal changes in the down- 
wash field due to the rolling -up process, and (3) to suggest a simple 
first approximation to the downwash at the tail due to a fuselage. The 
effect upon the downwash field due to substituting a lifting line for 
surface loading will also be investigated and an approximate method for 
taking this effect into account will be presented for wings of low aspect 
ratio . 


acn influence coefficients for a swept load line plus a swept trail- 

O IJL 

ing vortex sheet 

(These coefficients act as integration factors of the wing load- 
ing at station n to obtain downwash at position 

am^ influence coefficients^ similar to asn^ for only an unswept 

trailing vortex sheet (no bound vortex) 

h wing span measured perpendicular to the plane of symmetry, ft 

c local wing chord measured parallel to the plane of symmetry, ft 

g 

Cav average wing chord, 


PRINCIPAL NOTATION 


A 



mean aerodynamic chord 



c 


/ 1'^ 


local lift coefficient. 


local lift 
qc 


wing lift coefficient. 


lift 

qS 


lift -curve slope, per radian or per deg 


integration factor for interpolating downwash in the vertical 
direction 

empirical relation giving the effect of sweep upon the rolling-up, 
f = 1 -0.0075 (A°^ + 7°) 


strength factor of the tip vortices. 


^c 

r(n = 0) 


strength factor denoting loss of vorticity of the trailing sheet 
at span station n 

spanwise loading coefficient or dimensionless circulation along 

Q 

wing quarter-chord line, ^ 

G(*n) at span station t) = cos ^ 

o 

spanwise loading coefficient for 

op, 

— Gn where n refers to the 

Cl 

free -stream Mach numher 
free -stream dynamic pressure, Ih/sq ft 
radius of fuselage, ft 


or 


hV 


lit lift, ( ) 

\Ct cav/n 


or 


wing area, sq ft 

free -stream velocity, ft/sec 

downwash, positive downward, ft/sec 

right-hand Cartesian coordinate system with x positive down- 
stream and y positive to the starboard with the origin at 
the apex of the wing quarter-chord line (See fig. 1.) 


vertical distance in wing semispans measured from extended chord 
plane, positive upward 

inclination of wing from zero-lift attitude, deg 
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T 

cp 

n 

ftc 


av 

c 

f 


n,V 


s 

TE 


slope of camler line at trailing edge relative to free stream, deg 
-/l - 

circulation, ft^/sec 

angle of downwash, positive downward, radians or deg 

sweep angle of the wing quarter-chord line, positive for sweep- 
hack, deg 

tan-i 4eg 

tip chord 

taper ratio, 

root chord 

X y z 

dimensionless Cartesian coordinates, h^ 

longitudinal position at which sheet is essentially rolled-up 
into wing tip vortices 

yc 

lateral position of center of wing-tip vortex, — ^ 

h/c 

dimensionless longitudinal coordinate, measured from the lifting 
line (I - Ti tan a) 

trigonometric spanwise coordinate (cos"^ t]) , radians 
height above trailing sheet, ^ 
height above wing tip vortices, ^ 

Subscripts 

average 

tip vortices 

fuselage 

integers corresponding to span stations given by q = cos 

o 

or q = cos ^ (For n or v = 1, 2, 3> or It-; or = 0.9239^ 
0.7071, 0.3827, or 0.) 

pertaining to downwash at the sheet or displacement of the sheet 
wing trailing edge 
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PHYSICAL PROBLEM AND BASIC ASSUMPTIONS 


The physical picture is one of a lifting surface shedding a trail- 
ing sheet of vortices. As the trailing vortices are left farther behind 
the wing, the sheet of vortices is displaced downward in varying amounts 
depending upon the span station considered, that is, it assumes a curved 
shape. While this displacement is going on, the vorticity in the sheet 
is continually shifting from the sheet toward the tips or edges of the 
sheet. The lifting surface and the trailing vortex sheet are inclined 
with respect to the free-stream direction. 

Z 



The first assumption for the analysis will be that all of the chord- 
wise lift is concentrated at the chordwise center of pressure which will 
be taken as the wing quarter-chord line. Second, it will be assumed that 
the flow on the wing is not separated. Third, it will be assumed that the 
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downwash due to a symmetrical sheet can he approximated hy a horizontal 
flat sheet passing through the symmetrical sheet at the lateral station 
where the downwash is to he computed. It should he noted that at 



each T] station, the horizontal flat sheet is given a different vertical 
location and thus some allowance is made for the shape of the sheet. 
Fourth, it will he assumed that the vertical -longitudinal inclination of 
the system has no effect upon the downwash. Hence, the real system will 
he approximated hy a horizontal flat system passing through the real 
system at the downstream station, x, at which downwash is to he computed, 
as is shown helow. The coordinates of the real and substitute systems 



are shown in figure 1. It should he noted that these four assumptions 
are identical with those made hy Silverstein and Katzoff in references 1 
and 2. The first two assumptions are common in aerodynamics and the 
limitations are fairly well known for the higher aspect ratios. The 
first assumption will now he further investigated for wings of fairly 
low aspect ratio. 

Two wings having taper ratios of 0 and 1.0, aspect ratio equal to 
2.0, and sweep angle of 56° were investigated. Each wing was assigned 
hoth cotangent -type chord loading and uniform chord loading. The span- 
wise loadings were obtained from reference 5 » For each wing and chord- 
wise loading, the downwash in the wake, 6g, was computed with each of 
three alternative approximations; namely, the chordwise loading was 
replaced, respectively, hy a single lifting line, hy three lifting lines. 
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and "by five lifting lines. The strength and chordwise positions of the 
lifting lines were set by dividing the chord into equal segments and 
finding the lift and center of lift of each segment. Each lifting line 
was treated as a flat^ horizontal vortex system in the H = 0 plane. 

The downwash angles of the lifting-line systems were added and the sums 
are plotted in figure 2 for four spanwise stations. 

Figure 2 indicates that the single lifting line does not give 
accurate downwash predictions just aft of the trailing edge of the wing. 
The downwash fields for the wings of equal A are essentially the same 
one mean chord (one semispan) aft of the wing trailing edge. This con- 
currence at one mean wing chord aft agrees with the two-dimensional 
example of reference 1. For cotangent chord loading, the five-lifting- 
line method very nearly predicts c/a equal to unity at the wing trail- 
ing edge. This can be considered as a check to the approximation since 
the flat-plate downwash must be equal to a at the trailing edge. 
Examination of figure 2 shows that the curve of downwash obtained by 
using one lifting line is translated forward a nearly constant longitu- 
dinal distance from the curve of downwash obtained by using five lifting 
lines. In figure 2(a), this distance is one eighth the mean wing chord. 

In reference (l), contours of downwash angles due to a two- 
dimensional Clark Y airfoil section are compared to contours of down- 
wash angles computed for a lifting line at the c/4 point. If the lift- 
ing line is shifted back to the (3/8) c point, the shifted field agrees 
well with that of the airfoil section even very near the trailing edge. 
From this, it would appear that the downwash field due to surface loading 
might be well approximated for all wings by using a single lifting line 
with all longitudinal distances reduced by (i/ 8 )c 3 ^y^ or replacing t by 
T - (l/4) (c/b)g^y . It should be noted that this correction is of signifi- 
cance only in vicinity of the trailing edge. 

The third assumption has been considered by comparing the results 
obtained by using the assumption against results calculated for an ellip- 
tically shaped sheet whose ratio of minor to major axes was 0.4. At 
T] = 0, 0.383^ and 0.707^ the difference of the results was less than the 
differences found in the examination of the first assiunption. At 
T] = 0.924, use of the above third assumption did not compare well with 
the results for the elliptically shaped sheet. However, at low angles 
of attack, since the distortion of the sheet is small the downwash can 
still be computed at t] = 0.924. The fourth assimption has been checked 
by numerical computation for a 60° sweptback wing of aspect ratio equal 
to 3* 5* It was found that provided that e of the noninclined system 
is taken as W'/V rather than tan"^(w/v), the difference between the 
downwashes was less than the differences noted in examination of the 
first assumption. This appears to hold true up to about a = 20°. 

Thus, throughout this report, e will be taken as w/v and the subject 
is thus treated as if only small angles were involved. 
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It should also he noted that these four assumptions are commonly 
used in the calculation of wing span loading. As a result, the "non- 
rolling -up" system can he treated in a manner analogous to Multhopp* s 
(ref. 5) or Weissinger's (given in ref. 6) approximate integrations in 
their calculations of span loading. However, a principal problem not 
encountered in span-loading work is the downwash at arbitrary vertical 
locations. 

Generally, the amount of rolling-up present is so small that the 
foregoing assumptions are sufficient for good answers. However, as 
Cl/A increases, an increasing amount of rolling -up appears and a cor- 
rection must he made for this effect. The principal features of a 
trailing -vortex system where the rolling-up is conspicuous are, (l) the 
vorticity becomes vertically displaced and shifts outboard from the 
plane of symmetry, and (2) the wing tip vortices trail back approximately 
in a horizontal plane which is parallel to the free stream. The center 
of the sheet, however, is still displaced downward. As the vortex sheet 
is left farther behind the wing, the tips of the sheet roll up and form 
concentrated tip vortices. An outward motion of the vorticity in the 
sheet between the tip vortices results in less vorticity in the mid- 
semispan regions. These two changes in vorticity configuration can (in 
the main) be taken into account by making a fifth assumption, (l) a 
vertically displaced trailing flat sheet having a reduced amount of 
vorticity, and (2) a pair of tip vortices which "lie in a horizontal plane 
and whose strength is drawn from the sheet. With this arrangement, the 
sheet can be handled in much the same fashion as the flat sheet, that 
is, by using the first four assumptions. The tip vortices can then be 
handled as a separate computation. 



Sketch ( d) 

At various distances behind the wing, the rolling -up is in various 
stages of development. To obtain an accurate approximation, one should 

^ Wonrolling -up system assumes that the trailing vortex sheet has 
the same lateral distribution of vorticity at all distances behind the 
wing as at the wing trailing edge. However, it need not be flat although 
for determining downwash it is assumed flat. 
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consider the trailing system in longitudinal segments, each segment 
having a different amount of rolling -up. The downwash would then he the 
sum of the downwashes of all the segments. However, this involves an 
exorbitant amount of work and to obtain a practical solution, a sixth 
assimiption will be made. It will be assumed that the entire trailing 
system behind the wing is of one form, namely, the form which the real 
system has at the selected downstream location I . The substitute 
rolling-up system is then pictured as shown in sketch (e). 


C 



This sixth assiimption was examined by numerical computations for a 60° 
sweptback wing of A = 3*5 using the segment approximation. It was found 
that the results of the use of this ass\jmption were within the accuracy 
of the theory for this case. 

While the foregoing assmptions aid in simplifying the physical 
picture, additional information is necessary in order to calculate the 
effect of the rolling -up process. The relative strengths of the vortex 
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sheet as well as the tip vortex and also the position of the tip vortex 
for various distances behind, the wing must he obtained. From an analysis 
of the downwash behind a series of swept -wing plan forms obtained from 
large-scale wind-tunnel data, an empirical relationship was developed 
giving the approximate lateral position of the tip vortex. From this, a 
method is developed for obtaining the relative strengths of the tip vortex 
and flat sheet. 

As will be shown in the text, the use of such a simplified substi- 
tute system enables one to express the downwash due to a rolling-up 
system as being the flat-sheet results plus an additive correction 
(generally, fairly small) which, for the case checked, predicted the 
downwash well and goes to the right limit as the rolling-up becomes 
complete. 


ANALYSIS AND DEVELOPMENT OF METHOD 


The first part of the analysis is concerned with the flat-sheet 
proced\are, that is, evaluating the downwash using the first four assump- 
tions. The location of the wake relative to the tail will be considered 
and some assessment of the effects of the fuselage upon the downwash is 
to be made. In the second part of the analysis, the rolling -up of the 
trailing vortex sheet is considered. 


Flat-Sheet Procedure 


General calculation of downwash. - The downwash at a point (x,y,z) 
due to a swept-vortex system (assumed for the present at zs = O) is 
equal to the sirni of that due to the swept-load vortex (or bound vortex) 
and that due to the trailing vortex sheet. 
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The induced velocity due to an arbitrary elemental trailing vortex 
extending to infinity downstream (see sketch (f)) is given by Glauert 
(ref. 7). 


n P 

dw = •; (cos 0, + COS 0p) COS 

krth ^ ^ 


( 1 ) 


where © 2 — ^0. The vertical induced velocity at 
small element ds of the load vortex (see sketch 


(x,y,z) due to the 
(g)) is given by 


, Fh ds cos 
dw = a 


4jtr3 


( 2 ) 
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The total downwash due to the entire vortex system is equal to the 
integration over the wing span of the sum of equations (l) and (2) (which 
are converted to rectangular coordinates hy the relations indicated in 
the respective sketches). Further, the integral of equation (2) is inte- 
grated hy parts such that the new integrand contains the factor dF/dq. 
Then in terms of dimensionless relations, the total downwash can he 
written as: 


(q-q)G' (fi)dfi i 


Lsq G*(q) dfj 


w 

V 


1 


-1 


-1 


( 3 ) 


where for ^ > 0 ^ 


(Tl-ii) [t + ( |ii|-|fi|)tan A ] [ t + ( | t) | - T])tan A ] ^ -ti- |t) | tan^A - t tan ^ 


-\2 


■^S|-L 




[t + ( [T)[-'n)tan A] + n^/(cos^A) 


T1 - fj 


^[t + ( |T)|-|fj|)tan A]^+ + 


- \2 


+ (t| -fj) 


(M 


and for fj < 0, 


(Tl-n) [t + (UI-liiDtanA] [t+ (h|+Ti)tan A] f - ^ + [ti [tan^A + t tan A 

cos^^A 


_fi)‘ 


-•sia = 


[t + ( T 1 + Ti)tan A]^ + n^/(cos^A) 


-/[t + ( |r| (- [ri [)tan A]^ + (t) - fj)^ + 


XT+ (h|-Tl)tanA](T tan A+ h|tan^A + r^) ^ [r ^ ( h | + n)tan A](t tan A -h [n|tan^A-T i) 

^ + (hi -Ti)tan Af + flV(cos^A) [t + (h| + il)tan A] ^ ( cos^a) 

J~{r + |ti |tan A)^ + t|^ + 

^ - ii 

+ (ti - fi)2 
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Both integrals of equation ( 3 ) are, in general, difficult to inte- 
grate By analytical means. The numerical integration method used herein 
is that employed By Both Multhopp and Weissinger. In fact, equation ( 3 ) 
appears in reference 6 for the case of fi = 0. As was done in refer- 
ence 6 , the downwash integrals can Be written as a summation of products 
of mathematical coefficients, agn (avn in refs. 3 and 6 ), and the span- 
wise loading, Kn = (czc/Cl Cav)^- i^ shown in Appendix A, the n\ameri- 
cal integration of equation ( 3 ) can Be put into the form (using a four- 
term expression) . 


where for n = 1, 2, 3, and 4, the Kn’s correspond to values of 
czc/CLCav at T] = 0.924, 0.707, O. 383 , and 0, respectively. The agn's, 
like the Lg,^ in equation (3), are solely dependent, for a given sweep 
angle, upon the location of a point (T,r|,fl). Thus, after the agj^*s 
have Been evaluated, one may compute the downwash using any desired span 
loading. To facilitate computations, the agn's (from general equations 
of Appendix A) have Been computed for points at q = 0, O. 383 , 0.707, and 
0.924 lying Behind the quarter-chord line and for two vertical locations 
relative to the vortex sheet, (2 = 0 and (2 = ±0.5, that is, at the sheet 
and one -half semispan above or Below the sheet. The computed values have 
Been plotted in figures 3 and 4. Thus, given a plan form, span loading, 
and the desired longitudinal position, one obtains values of agn from 
figures 3 and 4 and applies equation (6) to obtain the downwash. 

It is obvious that the accuracy of a summation depends upon how many 
terms are considered. For the four-term s\immation used herein, accept- 
ably accurate answers are obtained in most cases without an unreasonable 
expenditure of labor. However, for cases in which the spanwise loading 
differs from that expressible By a four-term sine series, the four-term 
summation may not Be acceptably accurate as it would tend to gloss over 
such changes in span loading. The derivation of the agn’s in Appen- 
dix A has Been left in a fairly general form so as to allow the reader 
to compute the agn's for summations involving more points across the 
span. 

Choice of vertical coordinates .- Three possible vertical coordinates 
are the parameters (2, i, and Z. The vertical position of the downwash 
point for these three is measured respectively from (l) the trailing 
sheet, ( 2 ) a horizontal plane (parallel to free stream) through the apex 
of the load line, and (3) the extended chord plane. Each has some advan- 
tages. With (2, the downwash field is symmetric about the value (2=0 
(for a flat sheet), or about the trailing sheet. Also, computations are 
simplified with this parameter and the downwash varies linearly with C^. 
On the other hand, with the downwash field is referred to a fixed 
coordinate system, independent of angle of attack and sheet position. 

With Z, the downwash field is referred to the coordinate system (extended 
chord plane) of the airplane But is dependent on angle of attack. Thus, 
with Z, the downwash field is described relative to the tail plane. 



( 6 ) 


n=i 
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The relations between ^ , and Z are : 


« = 5 - 5g(Tl) 

Z = n+ i tan a 

o 

= ^ + 5 tan a 
= ^4- (t + T] tan a) tan a 


( 7 ) 


When the dovnwash has "been computed for various H’s and the dis- 
placements, have "been evaluated, the field can be plotted against 

a choice of Q, or Z. 


Lateral interpolation of downvash .- The ag^ values of figures 3 
and 4 allow a direct evaluation of dovnwash at four span stations. Some- 
times it is of value to know the dovnwash at other span positions or to 
plot a more accurate lateral variation of dovnwash. For these purposes 
a lateral interpolation formula is convenient. 

The product e( 9 ) sin cp can be expanded in a Fourier series, the 
coefficients of which can be numerically evaluated in terms of the four 
known values of Then (for symmetric distribution of dovnwash) 


e(cp)sin cp = 


sin l^j^cp 


,odd 


€n sin cpn 


sin 


n=i 


( 8 ) 


where cp = cos“^ Then the downwash at a given t] position can be 
expressed as the sum of products of tabulated numbers and values of the 
known downwash . Thus 


€ = Hid + ^ 2^2 + ^363 + 


( 9 ) 


where the H's are tabulated in table I for t] = O.098, 0 . 195 ; O.29O, 
0.556, 0.831, and 0.981, and are the known values of downwash at 

Ti = 0.924, 0.707, 0.383, and 0. 

Examination of table I shows that in the range of ti from 0 through 
0.556, Hi is very small, less than 4.5 percent of the sum of Hg, H3, 
and H^. For this range of t], one can simplify the calculations by 
letting €i = £2, then equation (9) reduces to, for 0 < q < 0 . 556 , 

e = (Hi + H2)€2+H3eg + (lO) 

This method of interpolation, in effect, puts a curve of the form 


e(T)) - eg + egTi^ + + ggT)® 
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(where the e’s are constants) through the known values of e at the 
four regular span stations. Of course, equation (lO) goes through only 
three of the four known values. 

The lateral interpolation formula applies for a given Q along an 
arbitrary ^(q) curve (e.g., at constant I or at a constant t). 

Vertical interpolation of downwash .- The asn values of figures 3 
and U allow a direct evaluation of downwash at fl = 0 and ±0.5. To 
present asn charts for many values is too cumbersome. However, 
with only two values of available, an interpolation procedure is a 
prime necessity. 

Reference 8 presents a simple method in which downwash, for small 0, 
is expressed as a Taylor series of |f2|. Reference 7 contains a simple 

equation of downwash as a function of ft for large ft values. With 
these functions of ft and the computed value of downwash at ft = ±0.5, 
a fitted function of ft is developed in Appendix B that approaches the 
correct functions at low and high values of ft and fairs through the 
computed value at ft = ±0.5* 

The vertical interpolation function is given as 

4 

e = C 3 ^€(t,t],o) + C 2 G^t,ti, |^± ^ ^ CnKn (H) 

n=i 


where C,, C 2 , and Dn are tabulated in table II for several ±ft values; 
€(t,t],o) and 6(t,t],±i/2) are the values of downwash computed by equa- 
tion (6) at ft = 0, and ±l/2, respectively. 

Vertical displacement of the sheet . - The vertical displacement of 
the sheet is given by the integration in the longitudinal direction (for 
constant q) of the downwash in the sheet from the wing trailing edge to 
the T and q position of the downwash point. Thus, for a given q, 

^s(”^) = ^TE ■ J' eg[Trji,^s(T,ji)] dTij, (12) 

Tte 

where vertical displacement of the sheet, Tiji is a dummy 

variahle of integration, and is the vertical displacement of the 

wing trailing edge. Now ^2 is defined as ^ then 

integrand corresponds to the value (2 = 0. With the flat-sheet assump- 
tions, €s (the downwash angle at the sheet, i.e., (2 = O) is independent 
of any vertical parameter and can he integrated to evaluate bs(T) . 
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The integration in equation (12) can he simplified considerably hy 
approximating eg with an integrahle function of Tip. An analytical 
study of figure 2 indicates that the downwash behind the surface -loaded 
wing varies closely as 1 /t^, also that the single -load-line method 
underestimates the downwash near the wing. The simplicity of the 1 /t^ 
behavior suggests a very convenient curve -fitting function for downwash 
that can in addition be made to correspond with the downwash near the 
wing as predicted by surface -loading methods. Thus it will be assumed 
that the downwash is given by the following function (fi = 0, and using 
the dummy variable Tpi) 


^s(’’’t) = + 

(Tip - 


( 13 ) 


where Ei, Eg, and j are undetermined coefficients. For the determina- 
tion of Ej^ and two conditions are given by 


at Tip - Tpg, 


^Te 


at Tp = T 


€s(Tt) = known downwash at the ? 

position at which dis- 
placement is to be 
computed 


( 14 ) 


With and E 2 determined it remains to evaluate j. Comparison of 
the results of several values of J with the downwash fields given in 
figure 2 shows that for j = c/2h very good agreement is obtained even 
with Gg taken at a distance of two semispans aft of the wing trailing 
edge. 


With the value j = c/2h and the determined values of and E^, 
equation (I 3 ) becomes 


€ 

s 


(Tp) 





a, 


TE 



( 15 ) 


Inserting equation (15) into equation (12) evaluates the vertical 
displacement. 
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S 




T 


“ "^TE 







oc. 


TE 


where Gg 


is computed at any t] 


position hut with (2 = 0. 


(16) 


For a wing having an airfoil section such that the load line is at 
the one -quarter chord and with equation (l6) simplifies to 




(t + T] tan a) a 



^ + T] tan A 


tan 



( 17 ) 


It should be noted that the last term of equation (l?) is negligibly small 
for many practical cases. 


Equation (l6) expressed as vertical displacement of the trailing 
sheet or wake from the extended-chord plane is given by (see eq. (7)) 



For a wing having an airfoil section such that the lifting line is 
at the one -quarter chord and for = a, equation (l8) simplifies to 
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( 



) 


(tan a - Gg ) + 




) 


c 


(tan a - a) (19) 



For many practical cases, the last term of equation (I9) can be 
neglected. 


Correction for the Effect of the Rolling-Up of the Sheet 


General analysis . - If the vorticity shed from the wing is now con- 
sidered to make up a weakened sheet and two tip vortices, as discussed 
earlier, then the correction, Ae, to be added to the downwash calculated 
in the foregoing sections can be written as 


where Cyg is the downwash due to a weakened vortex sheet; is the 
downwash due to a pair of tip vortices. The weakened sheet and the tip 
vortices extend from a longitudinal position corresponding to the 
quarter-chord point of wing tip, downstream to infinity. The quantity, 
G,pg, is the downwash due to the portion of the flat sheet aft of the 
quarter chord of the tip. Letting 


which gives Ae equal to the downwash due to the tip vortices minus the 
downwash due to the loss of vortex strength in the trailing vortex sheet. 

The downwash due to only a trailing sheet extending from the quarter 
chord of the rearmost wing section can be deduced from equation (l) by 
considering an unswept wing. Replacing i by ( I - tan A) then refers 
the coordinates to the quarter chord of the wing tip. Following the 
same mathematical procedures as in Appendix A, it is shown in Appendix C 
that the downwash due to such a trailing sheet can be reduced to 


■ ^WS ^TS 


€t - G.pg - €yg 


Ae becomes 


Ae = Gq - e,p 


(20) 
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( 21 ) 


n=i 


The difference between ^fiat sheet ^weakened sheets 

sidered as the loss of downwash due to the weakening of the sheet. This 
weakening can he taken into account hy multiplying Kn in equation (2l) 
hy a term called a strength factor. 

Thus, reduces to 


where Fgj^ will he called the strength factors of the sheet and denotes 
a loss of strength in the trailing vortex sheet. Thus, is the 

proportion of vorticity left in the sheet. The a^T^^’s are plotted as 
a function of (^ - tan A) /p in figure (5). 

The downwash due to a pair of tip vortices is equivalent to the 
downwash due to a trailing sheet behind a uniformly loaded wing of span 
T^gh. Thus, €(, can he derived from the preceding work. The circulation 
of the tip vortices is where Fc denotes the proportion of 

wing vorticity in the tip vortices. Expressions for a/f^Cl are 

derived in Appendix C and values are presented in figure 6. This param- 
eter is plotted as a function of (i/3ti„)(| - tan A) for various values 
of and ^cAc* quantity, (vertical height relative to the 

tip vortices), is equal to ^ - a tan A or + Cg - cc tan A, where 
-a tan A is the vertical location of the tip vortex which will he pre- 
sented shortly. 

As yet, Fgn and F^ have not been determined and, as was stated 
earlier, one must detennine the locations of the tip vortices before 
proceeding to find the strength factors. 

Location of tip vortices.- Reference ^ presents an approximate 
(curve -fitting) equation for the lateral position of the tip vortices. 
This equation is given hy 


4 



( 22 ) 


n=i 



( 23 ) 
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where 


I 


c 


5.05 (1 - Tic 

'^oo 

° ^ lim ~ 

Vl - cos cp 


(2U) 


where = 1 /k^ and represents the asymptotic position of and 

indicates the start of the rolling -up process. In the last part of 
Appendix C, is reduced to the more convenient fom 


5 


c 




0.7315 Ki - 0.3959 K2 + 0.3030 K3 - O.UOO 


(25) 


It should be noted that equations (23) and (24) were derived for a wing 
having an unswept trailing edge. 

It has been found from experiment that, for a swept wing, the inward 
movement of the tip vortices is much slower than is indicated by equa- 
tion (23). The reasons for this slowness are not clear. However, the 
lateral positions of the tip vortices appear to be strongly dependent 
upon the sweep of the wing trailing edge. 

The following two considered opinions of the actions of the air flow 
behind such a wing are given here. First, in the region between the two 
swept wing panels, the vortex sheet (principally near the plane of sym- 
metiy) is above both the load vortex of the wing and the wing tip vor- 
tices and thus is subjected to an inward velocity component. This inward 
velocity tends to keep the vorticity in the midspan region out of the tip 
vortices and thereby increases the roll-up distance. A second action 
concerns the wing tip. For a sweptback wing, an outward velocity over 
the wing tips is generated due to the lateral pressure gradient resulting 
from the staggering of wing sections. This flow over the wing tip is here 
assumed the primary action that results in a further outboard location of 
the tip vortex relative to that of an unswept wing. Similar reasoning 
leads to converse effects for sweptforward wings. The velocities involved 
in the above phenomena are difficult to determine. A problem remains, 
however, in that an expression for (compatible with experiment) must 
be found. 

The problem of theoretically determining r]^ for a wing with swept 
trailing edges is even more difficult than for an unswept wing. There- 
fore, an empirical equation will be used herein. One approximate method 
that takes into account the effect of the initial outward lateral 
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location of the tip vortices (due to the sweep) is to multiply (1 - 'Hcoo^ 
of equation (23) hy 


f + (l/lO) ( I - tan a) 
1 + (l/lO) ( i - tan a) 


where f is a function of the sweep of the wing trailing edge. At 
large values of I - tan A, this factor approaches unity and 
approaches I/K4. In the midspan region as mentioned previously, there 
is an inward velocity imposed upon the trailing vortex sheet which results 
in a longer roll-up distance. It is found that this can he taken into 
account hy multiplying the argument of tanh in equation (23) hy f 
where f here is the same as the previous f for convenience of empiri- 
cal evaluation. The effect of this multiplication is to increase the 
roll-up distance hy a factor of (l/f)®. Then, assuming the tip vortex 
to start at the quarter chord of the wing tip^ (Iq = tan A), equation (23) 
becomes 


f + (1/10) (1 - tan A) 

^ _ ) “hanVi 

^2 / e - tan A \2/3 

1 + (1/10) (1 - tan A) 

1 _L 1 UdilJ .1 

\ K4 y 




(26) 


where Iq " ^ given hy equation (25) for io ~ 

Equation (26) can he solved for f, since, for small values of the 
argument in tanh, the argument itself can he taken. A cuhic equation 
results from which f can he determined from experimental measurements 
of . Five wings with sweep angles of the trailing edge ranging from 
-Ul.3° to 51.9° were used in determining f. To the precision of the 
experiments, f is given hy, 

f = 1 - 0.0075 (Aijig + 7*^) (27) 

It has been asstuned that f is independent of aspec^ ratio. It 
should he noted that f does not quite reach unity at A*^g = 0. Thus 
equation (26) does not quite reduce to equation (23) for wings with 
unswept trailing edges. The difference, however, is small and is not 
considered important. 


^This is in contradistinction from reference 4 which assumes the tip 
vortex to start at the wing trailing edge. 
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It should he realized that other expressions for t],, can he 
obtained. Since only the numerical value of T|g is important, the use 
of other expressions of iic will not affect the material following. 

Later in the discussion section, it is shown that values of from 

equation (26) compare well with the limited available experimental 
results. 

The vertical position of the tip vortices can be determined by using 
the fifth assiimption in the previous section. Then, since the tip vor- 
tices are assiimed to start at the quarter chord of the wing tip, the 
vertical displacement of the two tip vortices is approximately given by 
the vertical position of the quarter chord of the wing tip, or 


= - a tan A 


( 28 ) 


It should be noted that far behind the wing, after the rolling-up 
process is essentially completed, the influence of one tip vortex on the 
other causes a displacement that varies linearly with I . The use of 
equation (28) therefore is restricted to locations near the wing. 

Strength factors of the trailing sheet and wing tip vortices . - The 
rolling-up trailing -vortex system is greatly simplified with the vortex 
system divided into two parts, (l) a pair of rolled-up wing tip vortices, 
and (2) a vortex sheet stretching laterally between the two wing tip vor- 
tices (see sketch (d)). The problem is to determine what proportion of 
the total vortex strength each should have. A method attributed to Lotz 
and Pabricius in reference k (given originally in ref. 9) is readily 
applicable for unswept wings. A modified and somewhat simplified pro- 
cedure of this method is developed here for swept wings. 

The basis of the method depends on two vortex laws applied to the 
rolling-up trailing-vortex system extending downstream from each half of 
the wing. These laws will apply aft of the quarter-chord point of the 
rearmost wing section where the trailing system may be considered free 
and two-dimensional. These laws are: 

(1) The total vortex strength shed from each half of the wing is 
invariant with distance downstream. 

(2) The total lateral moment of the trailing vortices shed from each 
half of the wing (or the total lift impulse) is invariant with distance 
downstream. 

Let Fs(ti) denote a proportional loss of strength in the trailing 
vortex sheet, then 1 - Fg(ii) is the proportion of strength remaining in 
the sheet. Also let Fq denote the proportion of strength in the trail- 
ing tip vortex. 
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The vortex laws yield two results. First, since the total vorticity 
is constant, the amount of vorticity in the tip vortices equals the total 
amount lost hy the sheet, or 


Fc = 



r T. / ^ 1 

Fg(q) 

/ hq 

K(o) J 


dT] = Fg(o) 


(29) 


Second, the moment gained hy the tip vortices equals that lost hy the 
sheet, or 



Equations (29) and (30) are not sufficient to determine and 

Fs(q). However, the form of Fg(ri) can he selected so as to represent 
the physical actions of the rolling-up sheet to a reasonable approxima- 
tion. In any rolling -up problem involving two main vortices, the vortex 
sheet in the outer span regions is acted on hy the tip vortices more than 
hy those in the inner regions. The outer vortex sheet rolls up at a 
faster rate than the inner. For these reasons the decrease of vorticity 
in the sheet, denoted hy Fg(Ti) should become larger as t) becomes 
larger. A simple expression for Fg(T]) which approximates these phenomena 
is 


Fg(n) = Fg(o) + ( 31 ) 

Fs(t]) at = cos 9n can he written as Fg^, thus 

Fsn = Fs 4 + A cos^cpn (32) 

where is evaluated hy equation (29) . The parameter. A, can he 

evaluated hy using equation (30) . Thus 
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The integral can he evaluated numerically hy use of equation (C23) and 
has the value 


^2 kOt)_ ^ 0.1283 — + 0.1388 — + 0.0531 

K(o) 


K. 


K. 


(34) 


Combining equations (32), (33) > and (34) results in 


^sn ~ ^c 


- l)cos' 


1 + 


(f) 


0.1283 + 0.1388 Kg + 0.0531 Kj 


(35) 


The strength factor of the wing tip vortices is assmed to be 
dependent upon the wing loading distribution and the lateral position of 
the tip vortices, qg. The latter gives a measure of the extent the sheet 
is rolled-up. Near the wing, the tip vortex would (at least) equal the 
wing circulation at q = T|c^ 4he center of the tip vortex. At most, the 
tip vortex would equal the wing circulation at the span station corre- 
sponding to the inner edge of the tip vortex. The inner edge, here, is 
roughly estimated as twice the distance from the tip to thus located 

at q = 2qc - 1. 

Now, in the derivation of the sheet strength factors, some sheet 
strength is left in these outer regions. Therefore, the strength of the 
tip vortices will not be taken as the larger of the above two values, 
but will be taken as the stun of the vortices inboard to q = (3qQ-l) /2 
which is midway between q^, and 2q^ - 1. Thus (again, near the wing)-. 



Far downstream, when the rolling-up is completed 

Fc = 1 (37) 

Between the above asymptotic values, assume Fq to vary as the 
product of wing loading at (3qc"l)/2 and a linear function of q^. 
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Then 


Fc 



K(o) 


+ kgTlc) 


(38) 


The two constants, ki and are evaluated by the two conditions given 
by equations (36) and (3?)* The resulting equation for Fg is 



K 


( 


2 


) 


is the value of 


^av 


at T] 


3ne-l 

2 




(39) 


K 



is the value of 


c^c 


"^av 


at T) 





and 


n 


^00 


1 

K(o) 


It should be recognized that other procedures that might determine 
the strength factors more accurately can be used in the present calcula- 
tions. However, it should be borne in mind that a fair amount of approxi- 
mation in the strength factors can be tolerated since a small percentage 
change in the factors results in an even smaller percentage change in the 
computed downwash. In a later discussion, it is shown for an example 
wing that experiment and the above theory compare well with regard to the 
vortex strength in the tip region. The expressions presented also have 
the merit of being computationally simple. 

Vertical location of the sheet during the rolling -up process .- When 
the sheet is rolling up, both the downwash induced at the sheet (fl = O) 
and the location of the sheet, ^g, will be different from the flat-sheet 
results. As an approximation for the rolling -up sheet, the downwash due 
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to the rolling-up sheet is computed at (2 = 0 using the value given 

hy the flat-sheet method. Substitution of this value of downwash in the 
displacement equation evaluates the ^g for the rolling-up sheet. This 
value of is used throughout the balance of the rolling-up computa- 

tion^ replacing also the ^g in the flat-sheet downwash; that is, 0. 
for the flat sheet is then. the same as Q for the rolling -up sheet. 

The rolling -up correction procedure, including the recomputation of 
^s, is illustrated in a computing form to be presented shortly. The 
change in location of the vortex sheet is generally quite small in the 
inboard region, becoming larger in the outboard region. 

Special loadings .- Span loadings which have maximum values at span 
stations other than at the plane of symmetry cause the equations for the 
strength factors to break down and in some cases to predict the lateral 
tip vortex location as being outboard of the wing tips. Such span load- 
ings would have complicated rolling-up characteristics since each change 
in sign of the slope of a loading distribution indicates the possibility 
of a rolled-up vortex eventually appearing. At the wing-tip region the 
loading gradient is very large, approaching infinity at the wing tip; 
hence, the rolling -up is more pronounced at the tips. The loading gradi- 
ent in the wing region between two maximum values of loading in general 
never becomes comparable to that at the tips, or even large. Hence, the 
rolling-up in this region will be very slow as compared to that at the 
wing tips. Thus it can be assumed for these special loadings that only 
the vortex sheet outboard of the maximum loading positions will roll up. 

Then, when determining a rolling-up correction, the loading distri- 
bution to be used will be that with a straight faired line connecting 
the two maximum values of loading. The ordinates of the entire curve 
are then proportionately reduced so that the area under the curve is 
equal to the original unit area, as is illustrated below. 



The area under the new faired loading is made a unit area by divid- 

-Ati 

ing the loading distribution by 1 + / AK dq. The lift coefficient 
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used for the rolling -up method increases to the value 


^^roll 


( 


1 + 


r 


AK dTi Cl 


) 


(ho) 


o 


where the subscript, roll, indicates the value of lift coefficient and 
loading used to compute the rolling-up correction. 


While a large part of the downwash at the tail is due to the wing, 
the presence of a fuselage will alter the downwash to a sufficient degree 
so that downwash due to the fuselage should be considered. For computa- 
tional purposes, downwash due to a fuselage at the tail of an airplane 
can be separated into two parts: that due to the fuselage at an angle 

of attack and that due to the fuselage at zero angle of attack. Two 
effects are present for the fuselage at an angle of attack. The span 
loading for the wing-body combination will be somewhat different than 
for the wing alone. This altered loading, together with the correspond- 
ing distribution of "image" vortices will affect the downwash in the 
region of the tail. If the fuselage diameter is not too large, the load- 
ing can often be approximated by the wing-alone span loading. However, 
in view of the pronounced effect of span loading upon the downwash, it 
would be preferable to use the span loading corresponding to a wing-plus- 
fuselage combination. The other effect of the fuselage at angle of 
attack is that resulting from the "crossflow" component normal to the 
axis of the fuselage in the region of the tail. The importance of these 
effects depends largely on the ratio of fuselage diameter to wing span. 

A further discussion is given in reference 10. At zero angle of attack, 
there is a disturbance of the flow field due to tapering of the fuselage. 
The influence of the fuselage on the downwash is further complicated by 
flow separation on the after portion of the fuselage, but this effect 
will not be considered. 

The relative importance of the above effects depends largely on the 
particular configuration being investigated and the spanwise region of 
interest. For the experimental data available for the present investiga- 
tion, a calculation of the influence of the fuselage at zero angle of 
attack appeared adequate. 

Since the fuselage diameters of most airplanes are small compared 
to the length, slender-body theory will be used to approximate the flow 
near the fuselage. Equation (lO) of reference 11 gives the complex 
potential in the cross plane (normal to the longitudinal axis of the 
fuselage) of a slender body at an angle of attack. For a body of 
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revolution at zero angle of attack^ the real part of this equation 
reduces to 


cPo = ao In r 


where a.Q = Rf (dR^ /ix) . The radial velocity is then given by 


dr V r r 

The downward vertical component gives the downwash. 


Wf 

T 


-zRf (dRf /dx) 


•zRf (dRf /dx) 

2 . P 

y + 


( 41 ) 


where y and z are meas'ured from the center line of the body. It should 
he noted that at the after part of the body dRf/dx is negative for a 
tapering body and that the equations are for bodies of revolution. 


Effects of Compressibility 


The effect of compressibility on downwash is taken into account by 
use of the Prandtl-Glauert rule. That is, the longitudinal coordinates 
are simply stretched by a factor 1/3 and c^ is replaced by 3c 

The compressible and incompressible parameters are listed below. 


Incompressible 




or 


A 

A 

5 or T 


n 


c ( chord) 


Compressible 

(replace the incompressible values by:) 


A 


3 c 2 , or 3Cj^ 

PA 

Ap = tan-i[(tan A)/ 3 ] 
1/3 or t /3 
4 

c /3 



J 
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The displacement values, or Zg, are not affected by 3 except 
insofar as eg is affected by compressibility. This eg can be calcu- 
lated by making the above substitutions. The t values in the ^g or Zg 
equations are then not replaced by x/p. The downwash for the case of 
sonic speeds is calculated by taking the limit as p — >0. The resulting 
simplified downwash equations are given in the rear portion of Appendix A. 
The values of downwash when p = 0 are not here represented as the true 
downwash but rather as a simple limit point that aids in fairing a Mach 
number curve of downwash. 


Computation Forms 


It is expedient to summarize the present calculations in a simple 
computation form. The forms for the computation of downwash angle due 
to a flat sheet, and the correction due to a rolling-up sheet, are pre- 
sented as follows: 

Flat-sheet procedure. - 


Raf- sheet downwash Verticol displacement 


77= T = 


Column no. 


1 
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5 
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a 

a SI 
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Vertical interpolation of flat -sheet downwosh 
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Rolling-up correction procedure . - 


Required information 


Trailing -sheet loss 


Clq= 

Cl= 

a ( radians) = 

i- 

tan A = _ 

equation (26) 'J7c=_ 

tan A) = 


For span loadings with 


Equation (39) Fc = 


Equation (35) 
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Tip-vortex contribution Total downwash 
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DISCUSSION 


This section evaluates the prediction of downvash due to wings 
alone, due to wing-fuselage combinations, and due to rolling -up of the 
vortex sheet by comparison with experiment. Also, for a pair of plan 
forms, downwash contours predicted by flat-sheet theory are compared 
with the flat-sheet theory corrected for rolling -up effects. 


Comparison With Experiment 


Wing alone.- Comparison of estimations from flat-sheet theory with 
some measured values of downwash from reference 12 for a swept-wing plan 
form is given in figure 7. For this wing, the computed rolling -up cor- 
rection was very small; hence, only the flat-sheet results are presented. 
The computed values make use of the calculated wing loading distribution 
obtained from reference 6 and also of the experimental loading distribu- 
tion,® which was somewhat different. The computed downwash distributions 
due to both loadings are presented and the two span loadings are shown 
in figure 8 . The general conclusion is that the downwash prediction at 
the plane of symmetry is critically dependent upon the local loading 
distribution. This is because the downwash contributions of the vor- 
ticity on either side of the plane of symmetry are additive for a sym- 
metrical span loading. At outboard stations, the downwash is not so 
dependent on the local loading since the vorticity to the inboard side 
results in an upwash which tends to cancel the downwash from the out- 
board side. It is noted that at the outboard stations the experimental 
and theoretical vorticity (or Iqading) distributions are more nearly 
similar. The effects of loading distribution are most prominent at the 
sheet. Figure 7 shows the experimental downwash from conto\ir plots and 
the downwash computed using the experimental span loading to be in good 
agreement. The experimental and computed locations of the wake center 

are also shovn 1^0 “be in good agreement > 

^This ”experimen‘tal** span loading was esilmaied from consideration 
of experimental results of numerous wings. — 
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Wing plus fuselage.- In figure 9 ^ comparison of the computed down- 
wash (wing -alone flat -sheet results plus the fuselage "taper” correction) 
and of experimental downwash is presented for a particular configuration. 
The wing is characterized hy an aspect ratio of 2.88, taper ratio of 
, 0.625, and the quarter-chord line was swept hack 50^- Th^ wing was set 

at 2^ incidence relative to the fuselage which had a length of 3-02 wing 
semispans and a maximum diameter of 0.297 wing semispans. The downwash 
is given at 1.239 wing semispans behind the quarter chord of the c. 

The fuselage taper (dRf/dx) was about -0.2. In figure 9^ it is shown 
that the downwash due to this wing-fuselage combination is predicted with 
I reasonable accuracy by the results of the flat-sheet method plus the down- 

wash due to the fuselage taper. At t] = 0, the computed and measured 
downwash above the fuselage are in good agreement at the lower angles of 
attack. Below the fuselage, the maximum difference is of the order of 1^. 
The poorer agreement at a = 13^ may be due (in part) to a change in span 
loading from that which existed at the lower angles of attack. 

At T| = 0.383^ the computed and measured downwash is in similar 
agreement. At a = 13^^ the possible change in span loading results in 
some discrepancy. However, as was pointed out in the comparisons of 
wing-alone downwash, this change in span loading does not result in as 
large a change in downwash at t) = O .383 as occurs at t] = 0. At 
T] = 0.383, the fuselage correction is quite small. 

In summary, in figure 9 it is shown that good predictions can be 
made by adding the downwash due to fuselage taper with the downwash due 
to wing alone . 

I 

I 

Rolling-up correction . - In figure 10 are presented measured values 
of downwash from reference 13 together with three methods of prediction 
for a wing with A = 60°, A = 3-5, A = 0.25, and Cj^ ^ 0.5. The three 
methods are: flat-sheet theory, rolling -up corrected flat-sheet theory, 

j and a completely rolled-up theory. The latter is simply the downwash 

due to the swept hound portion plus two concentrated tip vortices located 
at Tic = 0.86^4-. It is seen that the rolling -up method agrees well with 
experiment and that the agreement is best at the more outboard and at 
the more rearward positions. It is interesting to note that only the 
rolling -up correction method agrees well with experiment. For | equal 
to 2.71 and 3.43, neither the flat-sheet results nor the results for the 
completely rolled-up vortex yield maximum downwash angles within 10 per- 
cent of the experimental values. At | = 2.02 and t) = O. 383, the experi- 
mental downwash angles appear questionable because the rolling-up method 
gave virtually exact agreement with experiment at Cl = 0.25* 

In applying the rolling -up correction, it is important that the 
lateral position of the tip vortices be closely approximated. In fig- 
ure 11 is presented a comparison of measTired (refs. 12, 13, and li<-) and 
computed locations of the tip vortices. The computed locations of the 
tip vortices are based on empirical equation (26) . Figure 11 shows good 
agreement for a number of downstream positions for several swept wings. 
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The agreement is within the precision of the experiments. It should he 
noted that these are the same wings which were used in determining the 
empirical constant^ f, at one-half semispan behind the quarter chord of 
the wing tip. 

It is also important that the method yield a reasonable estimate of 
the strength factor of the tip vortices. It is difficult to determine 
the precise strength of the tip vortex from experiment because the vortex 
sheet is connected to the tip vortex. An approximate check can be 
obtained by determining the total vorticity in the region of the tip vor- 
tex. The experimental total can then be compared to the theoretical 
total. To check the totals for the 60° swept wing, the downwash contour 
plots of reference 13 were used. 

Very near the tip vortex (just outside of the maximum g values 
where the sheet contributes little) , the downwash due to the tip vortex 
is approximately given by 


Let 


w ~ *^vortex 
V it(Tic-n) 


^vortex 

(Cl/2A) 


Values of p ranging from 0.57 'to 0.63 were obtained from the expert- 
mental downwash contour plots for I = 3-i^3 and a = 12° by using various 
T]’ S . 


The theoretical total is taken as the theoretical plus the 

theoretical amount of vorticity left in the weakened flat sheet^in the 
region between t) = (3"nQ-l)/2 and 1.0. At i = 3*^3 and a - 12 , the 
computed Fc is 0.U8 and the computed total in the weakened sheet from 
■q = ( 3 t 1(, - l) /2 to 1.0 is 0.136. This theoretical total of O.616 is com- 
patible with the experimental total and is considered a reasonably good 
check on this phase of the method. 


Comparison of the Downwash Due to a Flat Sheet 
and That Due to a Rolling-Up Vortex Sheet 


The foregoing has indicated that the rolling-up correction method 
gives an accurate picture of the downwash fields behind swept wings. As 
is shown in reference h, Cl and A are important parameters in the 
rolling-up. Hence, it would be of interest to compare the flat-sheet 
results and the rolling -up results for a few combinations of Cl and A. 
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For this comparison the following wing -tail combinations were 
selected. 



A = 2.0 As 4.0 

A= 40* A= 40® 

X=0.5 X=0.5 

Sketch (i) 

Contoured downvash fields at the tail location as predicted by flat-sheet 
theory and by the rolling-up corrected flat-sheet theoiy are shown in 
figure 12 for two values of 0.4 and 0.6. 

In figure 12(a) (A = 2.0)^ the rolling -up is prominent and three 
principal effects are noteworthy: (l) in general^ there is an upward 

shift of the downwash field in the more outboard areas; (2) the magni- 
tude of the downwash around t) = 0 is reduced; (3) in the mid-semispan 
region (around t] = 0.5)^ the vertical distribution of downwash is more 
uniform than for the flat sheet. In general, the maximum values of the 
downwash for the two systems are not greatly different, but their loca- 
tions do differ appreciably. 

In figure 12(b) (A = 4.0), it is apparent that the amount of rolling- 
up present is quite small and the three effects mentioned above are 
scarcely discernible. In fact, these two cases. Cl = 0.4 and 0.6, could 
be considered as borderline cases. It is realized, of course, that near 
the tips (viz, t] = 0 . 9 ) ^ the rolling-up may have a sizable effect. The 
discussion here is limited to the more inboard locations as shown in the 
figure . 
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For other wings with the same Cl^ sweep, tail position, and taper 
ratio, calculations indicate that for aspect ratios about one (or less), 
the trailing -vortex system is, in effect, rolled-up into two tip vortices. 
For aspect ratios larger than four, the rolling -up present is even less 
than that present in figure 12(h) . 

Thus, the general range of aspect ratios for which the rolling -up 
correction should he computed is roughly hounded hy A = 1.0 and A = ^.0 
at Cl = 0.6; that is, (3/2 )<(a/Cl) <6. 


CONCLUDING REMARKS 


An influence -coefficient-type method is presented for the rapid 
estimation of the nonrolling -up downwash fields behind swept -wing air- 
planes. Using similar techniques, an additive correction for the effects 
of rolling -up is also presented. For the cases compared, the downwashes 
predicted hy the above procedures agreed well with experiment. 

To facilitate computations, charts and graphs of pertinent functions 
are presented together with tested computing forms. It is believed that 
the procedures set forth will require less time than procedures employ- 
ing horseshoe vortices or discrete vortices. To obtain a simple and 
rapid method, a number of approximations and assumptions were made. Each 
approximation and assumption was investigated by various means and the 
range of applicability is discussed. Some findings of the present 
research are as follows. 

By approximating the longitudinal variation of downwash behind 
surface -loaded wings by a simple function, a very simple expression has 
been derived for the vertical location of the wake center. It is shown 
that the location of the wake can then be written as a linear function 
of the downwash at the center of the wake . This downwash is easy to 
determine and thus the wake location can be detemined very rapidly. A 
comparison of experimental wake locations and computed locations indi- 
cates that satisfactory predictions are made. 

In the mid-wing region (around V[ = O) , it is found that the com- 
puted downwash near the wake is critically dependent upon the span load- 
ing used in the calculations. Thus, one should obtain the best available 
span loading before computing the downwash at the tail. 

The experimentally determined paths of the tip vortices trailing 
behind several wings have been considered. It has been found that wing 
sweep had an appreciable effect upon the mechanics of the rolling-up 
and slowed the inward motion of the tip vortices to a considerable 
extent. An empirical correction has been developed which allows one to 
determine the tip vortex locations with due allowance for the effect of 
wing sweep. 
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Aside from an effect of the fuselage upon the wing span loading, it 
appears that an important effect of the fuselage upon the downvash at 
the tail can he considered as an additive correction to the wing-alone 
downwash for wing -fuselage comhinations . This effect (due to the taper- 
ing of the rear portion of the fuselage) appears to he valid for comhina- 
tions wherein the diameter of the fuselage is fairly small compared to 
both the length of the fuselage and the span of the wing. The correction 
(obtained by slender-hody theory) is expressed in a simple form and has 
been shown to be in good agreement with experimental results on one air- 
plane model. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif . , 


for Aeronautics 
Sept. 16, 19 ^k 
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APPENDIX A 

INTEGRATION OF EQUATION (3) TO OBTAIN DOWNWASH 


A general equation for tlie flat-sheet downwash is presented in the 
section "ANALYSIS AND DEVELOPMENT OF METHOD." The purpose of the present 
section is to reduce the downwash integrals of equation ( 3) to the form 
of aquation (6) • It will he shown in the later part of this section that 
in the special case of sonic speeds, the downwash integrals simplify con- 
siderably. The downwash due to a completely rolled-up vortex system can 
be considered as a special case of the flat sheet; namely, that of a 
wing with rectangular span loading. Thus, the first portion of Appendix A 
is concerned with the general flat-sheet system, whereas the second por- 
tion considers specific cases. 


General Solution 


Representing the arbitrary loading distribution by a series and 
replacing the lateral integration variable by a trigonometric variable 
allows the first integral of equation (3) to be evaluated analytically. 
However, the second integral of equation (3) can only be evaluated numeri- 
cally and may be evaluated in the same manner as that of reference 6. 

The first portion of the following will be concerned with the analytical 
integration of the first integral of equation (3). 

It can be shown that G'(cp) can be represented as 


m m 


G’ (cp) = ^ ^ 


n=i 

Then, the first integral becomes 


(Al) 


in 


(cos cp - ■q)G’ (cp)dcp 
^ J + (cos cp - 

o n— ± 




vhere 


m , V 

^ ^ i}) 


(cos cp - q)cos Pj^cp dcp 
+ (cos 9 - q)^ 


Pi 


=1 


(A2) 


(A3) 
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Equation (A2) gives the integration to the same degree of accuracy as 
the vorticity distribution is given by the series of equation (Al) . It 
remains to evaluate the integral of equation (A3) . Define 


^ oit (cos 9 - T])cos 11 ^ 9 dcp 


‘•1 Jt 


+ (cos 9 - Tj)‘ 


(A4) 


With the relation for odd ni 


cos ^ 3^9 = 2 COS 9 


-1 


L Zl(Wi-Z)l 
1=1 


cos(|ii-2Z)9 


equation (a 4) can be written as a recursion formula. Thus 








(Jmi+ 1 - V^.l) - £ i!("" i), 

1=1 


where 


(A5) 


’^^l 


2 cos^^9 d9 


(cos cp - 


(A6) 


Now can be expressed in a recursion formula, by dividing the denomi- 

nator of the integrand into the numerator, then 




n^r 

dl-2 


, for |i^ = even 


, for = odd 


+ 2nj - (ti^+d^)Jj3 2 


(AT) 


where n^r are hinomial coefficients where 

n = ^ 1^-2 
^1-2 


and 


r = 
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The 

Let u = 
equation 


integrations represented hy Jq and are foimd as follows: 

tan — then cos <p = ^ . dcp = « With these substitutions, 

2 l+u^ 1+u^ 

(a 6) for Pi = 0 takes the form 


Jn = ^ 


^ 

o [p+(q) (u^)+(r) (u^)^] y (u^) 


y (u^)d(u^) 
pf(q)(u^)+(r)(u^)^ 


These definite integrals are tabulated in reference 15 and give the value 
for Jq as 


J q+27^ 

where 

p = + (l-T))^ 

q = 2(n^ + - l) 

r = + (l+q)^ 

J 

A similar procedure gives 


(a8) 


(1/Vp) - (1/^) 


(A9) 


where the p’s and q’s and r’s are the same as those in equation (a8) . 

With tabulated values of J|j,^ (from eqs. (A7) through (A9))^ the 
application of the recursion formula of equation (A5) evaluates the 
desired integral of equation (A^) . For symmetric loading, values of 
1 ^^ for odd M-i are given in the following equations: 
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Ii = 1 “ {a^+ Ti^)Jo + TlJj^ 


I3 = - (1 - + hn ^) - (t]^ + (-3 + Jq + l(-3 + - 12^2^)J^ 

I 5 = [1 -12t]^ + i6ti^ + (12 -96ii^)ft^ + i6n^] - 

(n^ + Q^) [5-20 ti^ + l6r)^ + (20 ~96ri^)a^ + l6n^] Jo + 

r\[9 - 20t)^ + 16t)^ + (60 - 160t\^)Q^ + SOft'^] Ji 


ly = [-1 + 2kr[^ - 80 t|'^ + 64t)® - (2k - U 80 ti^ + 960 ri'^)fl^ + 

(-80 + 96Or\^)0.^ - 64«®] - (T)^ + n^) [ -7 + 56 t]^ - 112t)4 + 

64 ti 6 - (56 -6T2ti^ + 960Ti‘^)n2 + (-112 + 960ii^)n^ - 64 ft®] Jq + 
ti [-7 + 56 ii^ - 112ti4 + 64 ti® - (168 - iiftOTi^ + i344Ti4)n2 + 
(-560 + 2240 T] 2 )fl 4 _ 448 ^ 6 ] 


For numerically evaluating for high it may he simpler to 

use equation (A5) directly; that is, to tabulate numerical values of 
from equation (AT)^ then with a numerical value of to tabulate suc- 

cessive I|j.i* 

With defined^ equation (A3) becomes 


The coefficients bsn found for arbitrary Q, and q. The value 

of the integrals of equation (A2) is the summation of the arbitrary 
loading^ Gn^ and the bsn values. The integral of equation (A2) gives 
the downwash due to a continuous trailing vortex sheet at an infinite 
distance aft of any wing and also twice the downwash at the load line 
for unswept wings. For the case of elliptic loading^ the downwash is 
given directly by 1^ alone, a result derived by other procedures in 
the past (e.g., see ref. 7) • For symmetric loading 
only odd equation (A2) becomes 



m 



(All) 


Pl=i 


m+i 



(A12) 
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where with \±2_ odd 


Bsn “ 2bsn^ 

n ^ 

m+1 

2 

= isn^ 

n = 

m+1 

2 


(A13) 


The numerical integration of the second integral of equation (3) is 
obtained from the appendix of reference 6 by substituting Ls|i given by 
equations (4) and (5) for L(T],fi) in reference 6. 

In summary, for symmetric loading, the dovnwash integral of equa- 
tion (3) becomes 


where 


€ 


W 

V 


m+i 

2 

^ 3'Sn^n 
n=i 


M-i 

2 

11=0 

in 

2 T J ni+1 

Bsn = ^ ^ ^ 

|i^=i,odd 


B 

s 


m+i 
^ 2 


^sn 

2 


for n 


m+1 

2 




and is given hy equations (AlO) and (A5) 


f‘n|i 


^n|i ^ 

= ^n[i ^ 
^n|i 


„ / m+i 

for n f ^ 

2 


m+1 

n . 


n ^ 


n = 


m+1 

2 

m+1 


p > 0 


P > 0 

p = 0 
p = 0 


2 


2 


m 


^ ^ cos 

M.l=i,odd 

, „ nrt , piit 

where cp^ - — , and cp|_^ = :^. Table III lists values of f^j^ for several values of 


m+1 


(’l“T|a)[T + (’l-iiia)'fcan A] T(T]-fj|^ + t sin A cos A) 


AL, 




cos^ A 


S|J. 


J [r + (Ti-jii^)tan A]^ + 


(n + ^|i)[t+ ('H- fj^) tan A] (t+2t) tan A) (t) - fi|^ + T sin A cos A - 2 t] cos^A) 

^ + (t + 2t) tan A)^ cos ^A 

V [t + ('n-Tii^)tan A]^ + (il+Ti|^)^ + 


t(t) + t sin a cos a) (t + 2t) tan A)(r| +t sin A cos A -2q cos^A) 


+ t^cos^A 


+(t + 2t| tanA)^cos^A 


T - Ti. 


/i 


T + T) tan A)^+ + 




where 


V 


jin 

COS = cos CD,, 

M+l 


M and m. 




-P:- 

U) 
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A convenient closed form of fn^ is given by; 


[sin cp^ sin (m+l)cP|a(2 + cos 2cpn + cos 2<P|j,) - 

(m+l) (cos 2cpj^ - cos 29,^.)^ 

(m + l) cos cos (m + 1) 9|j. (cos 29^^ - cos 29^)1 


for 


= ^n 


2 sin 29ji 


for 


M = m 


f 


n|i 


n+u 

2(-l) sin 9n ^os 9|i 
cos 29n - cos 29^ 


Special Cases 


Compressib ility considerations and downwash at son ic speeds. The 
effects of compressibility, subject to the limitations of the linearized 
compressible flow equation, can be included in the previous work by the 

substitution of the parameters | and A^= tan'^^ ^ , for I and A, 

respectively. It can be seen that since the Bsn coefficients are inde- 
pendent of I and A, they are unaffected by compressibility and that only 
the Lom function is affected by compressibility. 

The limit value of the downwash at sonic speeds can be found by sub- 
stituting into equations (3), (^). and (5) the parameters t/P and Aq 
which replace t and A, respectively, then determining the limit as P^O 
With the limit 


lim 


cos Ap 


tan A 


(A15) 


<\< 


equation (3) at sonic speeds becomes 


1_ 

2it 


-1 


1 + 


T + ( I 'H I - I n I ) "tan A 

It + (h|-|ii|)tan A| 

[ T + ( I q I - T| ) tan A] tan A 
[t + (|ti| -T])tan A]^+n^ tan^A '■'o 

[t + (Iql + T] ) tan A] tan A 
[t + ( |q I + q)tan A]^ + Q^ tan^A J_^ 


T) - Tl 


G'(ri) dfj - 


T + 1 q I tan A T + (|Ti|-|fi|) tan A 


|t + I Til tan A| |t + ( |ti| - |fi| )tan a| 

T + I Tl I tan A T + (|Ti|-|fi|) tan A 


I T + I Tl I tan a| |t + (jri|-|Ti|) tan A 


G'(ti) dfi + 


G' (fi) dfj 


(Al6) 


The form of equation (Al6) can he simplified hy considering three longitudinal regions of 
downwash. These regions can he pictured as follows: 


o 






Region I 



-P" 

\J1 
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For the three regions, equation (Al 6 ) gives 
Region I: 




Region II (for sweephack) : 


w 

V 


2n 



T+ 1 rj 1 tan A 

-/ 

tan A 
r+ |t] I tan A 


tan A 


4-4 


+ (n-q)2 


G’(4)<i4 + 


^ [ T + ( 1 4 1 - 4 ) tan A] tan A G I — ^ ^ j 

\ ' tan A 1/ 

[t + (hi- 4 ) tan A]^ + (2^ tan^A 


2[r + ( I r) l+T)) tan A ]tan A G^- 


T + I Ti I tan A 


tan A 


[t + ( 14 | + 4 ) tan A]^ + tan^A 


» (AI7) 

; 
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Region II (for sweepforward) : 


T+ 1 T] [ tan A 


w 

V 


2jt 


tan A 


-1 


^ ^ 




G'(fi)dTi 


I T+ I T) I tan A 
tan A 


^ - T] 

+ (n-Ti)' 


G'(ri)dfi + 


2 [t + ( 


In I - Ti)tan A]tan A Q ^ 

[■’■ + ( 1^1 " n)tan A]^ + tan^ A 


Region III: 


2 [ T + ( I T) I + n ) A] tan A G I '’^'^ 1 ^1 ^ I 

\ I tan A I / 

[t+ (|ti| +n)tan A]^ + tan^A 


Tj - T1 


w ^ 1 


G* (n)dTi 


(A18) 


The symbol G( ± 
span station. 


T + 1 n I tan A j\ 
tan A ' / 


denotes the value of G at that 
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For synmie'tric loading. Region II given by equation (A16) simplifies 
(since G(-ti) = G(ti)) to: 

Region II (for sweepback) : 


w 

V 




T + q tan A 


1 

— < 

r 

tan A 

^ ~ ^ G’ (Tj)dTi + 

7t 

j 

T+ q tan A 

G^ + (Ti-n)^ 



tan A 



tan A 


r + 2r\ tail A 


_ tan^A (t + 2t^ tan A)^ + tan^A . 


T + q tan A 


tan A 


(A19) 


Region II (for sweepforward) 


V 

V 


1 
— 

7t 


t+t] tan A 


tan A 


r\ - r\ 


-1 




G’(fi)dfi + 


q - T] 


T+T] tan A 


tan A 


G^ + (n-^i)^ 


G’(ri)dTi + 


tan A 


T 4- 2t] tan A 


-T^ + tan^A (t 4- 2 j ] tan A)^ + ^ tan^A - 


T 4* T] tan A 
tan A 


(A20) 
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For the Region II for the case of swepthack wings, the downwash can he 
written as (symmetric loading) 


w 

V 


III 


1 

TC 


T+Ti tan A 


tan A 


-1 


a ^ ~ G’ (n)dfj 

+ (n-n) 


1 

Tt 


t+t] tan A 


tan A 




G’ (fi)dfj + 


tan A 

It 


T + 2t] tan A 


tan^A (t + 2t] tan A)^ + tan^A 


a^ 

T + Tj tan A 



tan A 



(A2I) 


The two integrals can he approximated hy replacing the loading distribu- 
tion near the wing tip hy a single vortex having strength given hy 
T T] tan A 


and laterally located so that the downwash at r) = ^2 = 0 

tan A 

due to elliptic loading at the wing tip equals that of a single vortex. 
Then the lateral location of a single vortex is at 


1 - 


T = 


T+T] tan A 
tan A 


cos 


tan A 
tan A 


ifi. 


T+T) tan A 


tan A 


(A22) 
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Then the downwash becomes 


w 

V 


III 


/ 


fl. 

t+t] tan A 


tan A 


" 


jt 


\/W* 

T+T] tan A 

V.l 



Lv 2 V 

tan A 

J J 


-i2 


+ 



^1. 

T+T] tan A 


tan A 

1 1 


T+T] tan A 
tan A 

)*r 


G 


T +11 tan A 


tan A 


tan A 


+ - 


T + 2t] tan A 


^ Lt^ + tan^A (r + 2 t) tan a)^ + tan^A_ 


T+T) tan A 


tan A 

(A23) 

The downwash in Region III can be evaluated from equation (A12), 


then 


For symmetric loading 


m 


III 




(A2k) 


n=i 


m+i 

2 


(0... ■ I 

n=i 


2BsnGn 


(A25) 


where the value of 2Bsn is the value of asn s-i t = oo given hy 
equation (Al4) . 

Downwash due to rectangular span loading .- For rectangular span load- 
ing the equations for determining downwash simplify considerably- Since 
the loading is constant across the wing span^ then in equations (3)^ 

and (5) one can substitute ”/^G*(fi)dfi = G(o)^ and plus and minus values 

o 


of unity for fi. The downwash equation becomes 


w 


/V 


2G(o) 


1 - X] 


(l-T))[T-(l-ri)tan A] 


T tan A - 


- T)) 
cos^A 


-f flVcOS^A 


fl + (l-Ti) +(i+T|) y [t - (l-Ti)tan A]^ + (l-T))^+ft^ 


(l+q) [t -( i-T])tan A] ^ (T+2q tan A)[t tan A - 2t| - (1-t|) /cos ^A] 
+ (1+^)^ (t + 2q tan A)^ + fl^/cos^A 


y [t - (l-q)tan A]^ + {l+T])^+n^ 


t(t tan A + ti/cos^a) , {r + 2i\ tan A)(x tan A -2ti + g/cos^A) 

T + fl^/cos^^A (t + 2t) tan A)^ + G^/cos^A 

y(r + q tan A)^ + t)^ + 

For n = 0 equation (A 26 ) simplifies to 

v/v 


2G(o) 


(A26) 


JL + _1_ + y[T- (l-q)tan A]^ + (l-g)^ y [t - (l-q)tan A]^ + (l+q)^ 
1-q 1+q _/■. ^ : : - + 


T(l-q) 


(t + 2q tan A) (1 + q) 


2 tan Ay (T-Hq tan A )^ + t|^ 
t(t + 2r\ tan A) 


(A27) 


VJI 
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CompressilDility is taken into account in the same manner as in the 
previous section. At sonic speeds for the same regions as before^ equa- 
tion (A26) reduces to 

Region I: 


v/V 

2G(o) 


= 0 


Region II (for sweephack) : 


w 


/v tan A 


2G(o) 2n 


T + 2t) tan A 


+ a^tan^A (t + 2ti tan A)^ + G^tan^A 


(A28) 


Region II (for sweepforward) ; 


w/v 

2G(o) 


if 1-3 1 + ^ 

+ (1-T))^ + (1+T))^ 


+ tan A 


T 

+ n^tan^A 


+ 


T + 2ri tan A 


(t + 2t] tan A)^ + _ 


Region III: 


w/V ^ J_ r 1 - T] ^ 1 + T] ~ 

2G(o) + (l+T])^ _ 


(A29) 


(A30) 


If the trailing vortices due to rectangular loading are not at the 
wing tip (t] = l) laterally located at then the dovnvash is 

obtained hy substituting ^ ^*or G(o), t], and 

'He *Hc ^c 'ic 

respectively, in the previous equations- 
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APPENDIX B 

VERTICAL INTERPOLATION OF DOWNWASH 


An approximate vertical interpolation of downwash can "be found by 
using the method of reference 8 for small values of Q. In the present 
report, the interpolation formula for other values of will be based 
on a fitted function of fl that fairs through the two known values of 
downwash (ft = 0 and±0.5)^ approaches the function given by reference 8 
for small values of ft, and approaches the correct function for high 
values of ft (e.g., ref. Y). 

In reference 8 the downwash function is expanded in a Taylor series 
starting from the vortex sheet: 


e(ft) - e(0) + IftI + I ft^ 


(Bl) 


where the subscript^ indicates evaluation at H = 0. Assuming the 
trailing sheet extends fore and aft to infinity, the following relations 
are derived in reference 8: 


8ft ~ 


8 ^ 6 ^ _ ^ 
8ft^ ^ 


(B2) 


Now, for large values of ft, if t is small relative to ft, the 
downwash for any loading distribution is given by 


e(ft) 


£l 

2(ft^ + r\^)^ itA 


(B3) 


As ft becomes very large (compared to t]), equation (B3) becomes 


€(n) 

2ft^ itA 


(B4) 


The derivative of equation (BU) gives 


dft 


.A. ?A 

nA 


(B5) 


It is desired to curve -fit a function of ft that approaches equa- 
tion (Bl) for small ft, gives the known value of e at ft = ±0.5, and 
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approaches equations (B3), (B4) , and (B5) for large Q. Assirne the 
following function: 



using eq. (AlO)) to he proportional to this parameter for elliptic load- 
ing at T) = 0. 


The coefficients and C2* are evaluated by taking the first 

and second derivatives with respect to Q of equation (B6) , which, at 
ft = 0, give the values of equation (B2) , that is. 


d^G 

dq^' 


1 ^ 
2 dq^ 


(B7) 


The four remaining coefficients of equation (b 6) can be found by 
using four conditions which the equation must satisfy. One condition is 
that it pass through the known value of e at ft = 1/2. Another is that 
it pass through the value of e at ft = 2.0 given by equation (B3).^ A 
third condition is that at very large values of ft, the downwash must be 
the value given by equation (B4) . A fourth condition is that the slope 
of the curve at very large ft be that given by equation (B5) . The four 
coefficients (A^^, A2, A3, and A4) so determined will be in terms of €q, 


-i 


d^e 


d^ 

dq^^ dq' 


and — . The next step is to evaluate the two deriva- 
jtA 


Ct 

tives and — . 

A 


To find 


d^G . 

dqS- 


Now 


d^G(cp) dG(cp) 




dcp 


d^T d^G(9) 
dT)^ dCp^ 


dq ) 


(b8 ) 


It should be noted that at two semi spans from the sheet the down- 
wash is essentially independent of span loading (e.g., ref. 7 , p. I65); 
however, the downwash still has enough magnitude to make it useful in the 
curve fitting. 
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where 


Cp = COS“^ T] 


Then 


d9 

1 

d^9 

dq 

” ^ 

sin 9 



cos cp 
sin®<p 


d^G(cp) ^ _ cos cp dG(<P) 1 d^(y) 

dT)^ sin®cp dcp sin^cp dcp^ 


(B9) 


(BIO) 


The derivatives of G(cp) are obtained from equation (Al) . Then equa- 
tion (B8) becomes 


d^G(cp) 

dq^ 


m 


-I 


Vr n 


n=i 


where 


m . 

2 y f cos [1(f) ^ \ n sin u cp 

= — 7 ) ( ^ 1^1 sin | 1 , cp^ )— i LJl 

\ tan cpy J sin^cp^ 


. n-v 


(-1) sin <p^[2 sin^cpv+ (cos cp^ - cos cpy)cos 9 ^] 


=( 


sin3cp^(cos cpn - cos 9^)^ 
m(m+2)sin^ 9v -3 cos^ 9^ 


3 sin‘^9. 


> n ^ V 


n = V 


(Bll) 


> (B12) 


For symmetric loading. 


d^(9) 

dq^ 


m+i 

2 

n=i 


(B13) 
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where 




n 




m+l 

2 



n 


m+l 

2 


For m = tbe 


Zv values are given as follows: 


(Bll^) 



0 

0.383 

0.707 

0.924 

1 

-1.7934 

3.5149 

-26.5029 

100.7692 

2 

5.6568 

-16.5757 

40.0000 

-19.1087 

3 

-25.2346 

27.2308 

-19.6374 

-20.4839 

4 

21.0000 

-13.5139 

2.8284 

11.7199 


With d^c/dTi^ evaluated, the next step is to evaluate d^e/dt]^. 
The downwash of the two-dimensional trailing vortex sheet is given hy 
equation (A2) . For n=0, equations (A2) and (A3) simplify and become 
(e.g., ref. 6 , eq. (AI9)) 


e(9) 



sin sin 

sin cp 


n=i M-i”! 


Then similar to equation (BIO) 




cos cp de 1 d^e 

sin^cp dcp sin^cp dcp^ 


(BI5) 


(BI6) 


The derivatives occurring in equation (BI6) are obtained by differ- 
entiating equation (BI5); then after summing up several resulting 
summations, equation (B16) becomes 


d^e^ 

dT)^ 


m 


I 




n 


(B17) 


n=i 
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where 


^ — [15 - (m^ + 2m + 12) sin^cp^] 

n=V 4 sin^cpy 


Evn - 


-3 sin cp, 


n 


n^v (m+l) sin^ <P-y(cos Cf^-cos 9^) 


(m+l)^(-l)'^ ^ cos cpy 

sin^cpv 


(Bl8) 


(m+l)^(-l)^“'' 2[l-(-l)^-"']sin^q), 

+ 


cos - cos cp^ (cos cpn ~ 


For syminetric loading 


where 


d^e. 


dTi' 


m+1 

2 


Evn'^n 


n=i 


J 


(B19) 




^Vn= < 




n ^ 


n = 


m +1 

2 

m +1 


(B20) 


For m = 7^ the values are given as follows 


X 

0 

0.383 

0.707 

0.924 

1 

19.945 

-48.262 

241.673 

1174.743 

2 

-67.883 

128.450 

-237.586 

-1791.706 

3 

173.641 

-174.776 

81.739 

1890.134 

4 

-120.000 

89.112 

0 

-929.163 


The lift coefficient can also he given as a svunmation of 
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Thus 


Cl 

A 




(B21) 


For m = 7, the factors of Gn equation (B21) are as follows: 


n 

1 

2 

3 

4 

factor 

0.3006 

0.555^ 

0.7256 

0.3927 


As mentioned "before, the coefficients of equation (B6) (A^, A^, Aq, 
and A4) are in terms of e(0), e(l/2), d^o/dq^, d^e/dq^, and Cl/hA. It 
has just been seen that d^o/dq^, d^e/dq^, and ClAa can he expressed 
in terms of summations with Gn* Thus, equation (B6) can he written in 
terms of e(0), e(l/2), and summations of Gn* The equation 

can now he algebraically rewritten into the form 


e(fl) 


m+i 

= e(0) + C2 + £ 

n=i 


hnGn 


(B22) 


where the constants C^, C2, and Dn contain the J terms and 
the integration coefficients of d^G/dr)^^ d^e/dr]^, and Cl/a. Letting 

X = ^/G^/l+^^^, the constants are: 


C. = 1 


172+111 ^3 + 539+282 ^/5' ^5 


610+231/5 ^5 235+60 /F ^6 


8 


(B23) 


C 2 = ^(5 + 3/5)X®+^(35 + 17a/5)x^+25 (5+ 2/5)X^- ^ (15 + 5/5)X 


8 
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47 + 12 ^/5 v3 , 98+39 75+42 ^/T ^5 20+19 X® 

1 „2^ WsVT x= _ x" + x' - f X®') Ev + 

'2 4 4 ^ A/n 

'l±i^x® - 59i£L/ix-‘ + (20+9^5) x= - 3?ti5V 3 x^ 1 ^ 

_ 2 4 + J 

■- ?5 (20+9^) X® + ^ (89*38 v^) X‘ - f2 (110+49/5) X® H 
16 18 -LD 


(9+4-/5) X® 

16 


(4-n^) 

(4+n^)^ 


} 


( 2 sin cp^ 
m+1 ^ 



For m = 7, values of C^, Cg, and 
various values of H . 



(B25) 


are presented in table II for 




6o 


NACA TN 3346 


APPENDIX C 

DEVELOPMENT OF EQUATIONS PERTAINING TO 
THE ROLLING -UP CORRECTION 


This appendix includes the downwash due to a flat trailing sheet 
alone, downwash due to a pair of trailing wing-tip vortices, and a 
reduction to more convenient form of the "rolling-up distance" of refer- 
ence 4. 


Downwash Due to a Flat Trailing Vortex Sheet Extending 
From the Quarter Chord of the Wing Tip to Infinity 


The downwash can be obtained from eq\iation (l) by setting 


cos 0, = 


X - (b/2) tan A 


y (x - I tan a') + (y-y)^ + (z-zs)‘ 


(Cl) 


In dimensionless coordinates^ the downwash becomes 


(Ti-fi)G' (ii)dfi 


-'-I + (Tl-Ti)^ 


- F - 

2n J 

-1 


T1 - T] 


+ (ti-r) 


s\2 


^ - tan A 


•/ (5 -tan A)^ + (T]-fi)^ + 


G‘(Ti)dfi (C2) 


The first integral is integrated analytically in Appendix A, and is given 
by equation (A12) . 

The nimerical integration of the second integral can be performed 
by following the procedure of Appendix A. Let 


Lth = 





I - tan A 

(^-tan A)^ + (rj-Ti)^ + 


(C3) 
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then the second integral of equation (C2) 


becomes 


6" 


= cos CPjj^ = 


cos 


m+1 ) 


2jr 



(C4) 


where 


M-i 

2 

' 2(M+1) X ■ ^T,M+i-q) 

q=0 


(C5) 


where fnq values are given under equation (Al4) . 

For the case of q = ft = 0, the second integral of equation (C2) 
can be handled in much the same fashion as equation (A2) and with the 
same limitations. The vorticity distribution is given by equation (Al) 
which for symmetrical loadings becomes 


dG(cp) 

d(p 


m+1 


Z 

2 Ju .=: 


Hi(-l) 




cos 


UlCp 


|i 1 = 1 , odd 


m-i 



n=i 



|i^=i,odd 


sin cos ii^cp 


(C6) 


With equation (C6 ) , the second integral can be written as 



m-i 


I y 

L -1 (m+i)7r 

n=i _ un = 


sin [i^cpn / cos ii^cp dcp 


(C7) 


odd 
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Now, for ft = T] = 0 
/ cos iij^cp d<? = 


Jt 

r 1 

^ - tan A 

1 

cos cp 

{i -tan A)^ + cos^cp 


-(i)- 


= jt sin 1 ^, ( ^ ^ - tan a) 


cos |j.^cp dcp 


cos jij^cp dcp 


cos 


9a/U -tan a)^ + cos^cp 
(C8) 

The integral remaining in equation (C8) can he written as^ 


cos dcp 


Jt/2 


cos cp»/ (I -tan A)^ + cos^cp 


= 2k 


cos dcp 
cos cp^/l-k^sin^cp 


where 


k = 


// 1 + (I -tan A)^ 

This integral is evaluated hy a recursion method. 
For odd values of 

Hl-i 


cos n^cp 1 ^ -1 ^ -1 

= 2 cos cp - 

cos cp 


1=1 


cos (ni- 2 l)cp 
cosCp 


Hence j 


Hi-i 

2 


= 2% R/ 


1^1 • 


L U (hi-Z)! 
^ ^ ^ 1=1 


(C9) 


(CIO) 


(Cll) 


iThe integration of equation (C9) was obtained with the aid of 
reference l6 . 
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where 


- 2k 


cos ii^cp dcp 


cos cp y l-k^sin^<p 


(C12) 




cos^^ ^cp d<p 
y 1-k^sin^cp 


(C13) 


For = 1 and = 3^ equation (C13) integrates directly to 


Bo=- = K 

Ri = i [E - (l-k2)K] (ci4) 

k 


where K and E are complete elliptic integrals of the first and second 
kind, respectively. Higher terms of can now be obtained from the 

recursion formula: 


X^) ° (^x-2)k= 


W) 


(C15) 


With values of determined, the values of are obtained 

from the rec\irsion formula of equation (Cll). Then 


= 2kK 


Da = - - k^)K-4E] 

Ds = -%[(32 - 20k^ + 3k^)K - 4(8 - k2)E] 


yci6) 


3k 


D., = 


^[(512 - 608k^ + 2l6k^ - 15k®)K - (512 - 352k2 + 72k^)E] 


15k ^ 

where k follows equation (C9). 
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Thus the second integral is evaluated exactly for ft = ti = 0 and, 
in simmary, is given by 


1_ 

2n 


J (cp)dcp = 


m 




(m+l) 




4i(-i) 


(V) 


('v£) 

It(-l)V 2 ) _ 




> + 


|j^=i^odd 


m-i 


l—i (m+i)it Z_i 


HiSin nfpn 


n=i 


bl=i,odd 


rt sin - (^-tan A)D^^ 


> (C17) 


vhere the values are given hy equation (Cl6) and for higher values 
of \± 2 _ hy the recursion formula of equation (Cll) . 

Similar to the work of Appendix A the downwash, given hy equa- 
tion (C3), can he expressed as the sum of the products of influence 
coefficients and values of the loading distribution, or 


where 


m+i 


w 


a^n^n 


n=i 


M-i 


arpn - 


2Bsn 


-1 


2(M+i) 


fn|jAIjji|_i 


|j=o 


(Cl8) 


The coefficient, Bsn^ follows equation (A14) . 


Similarly, fj^^^ also 


NACA TN 3346 


65 


follows equation (Al4) 


ALrii = 


•n - Tlu 


i - tan A 


U-tan A)^ + (rj-fiJ^ + 


3 + n. 


1 - 


i - tan A 


1_ y(|-tan A)2 + (q+f )^)2 + 


T]j, = COS 


m+1 


For T] = Q = Oy the analytical integration which evaluates is 


2(1 -tan A) 


^Tn “ ®sn < 


S-tan A 

(m+l) 


m 

V T^ • / 

/ Hj^cp, n ^ 

^ |ii=odd 
m 

- ^ ^ , n = 


Hj^=odd 

where values are given in equation (Cl6). 


m+l 

~2~ 


(C19) 


Downwash Due to Two Trailing Vortices 


Let the lateral position of the two vortices he at fj = where 

is the fraction of wing semispan from mid-wing to the center of the 
wing tip vortex. The downwash equation is obtained by substituting 







/ G'(fi)dfi = G(0) 

p'lus and minus values of t)« 

0 

for T], and T)^^h for h in equat 

Then 



w _ G(o) 


i - tan A 

1 + 

V ~ 2n ' 


_ y (l-tan A)^ + (q^-'n)^ + 


T).+ Ti 


+ (q^+q)2 


1 + 


i - t an A 


y (l-tan A)^ + (nc+'n)^ + 


(C20) 
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Now, can be absorbed into coordinate parameters, tben 



(C2l) 


Roll-Up Distance, Iq " ^o 


The longitudinal distance at which the vortex sheet is essentially 
rolled-up is a function of wing loading, lift coefficient, aspect ratio, 
and sweep. The roll-up distance given in reference U is given by (when 
the coordinates are changed to the present notation) , 


I 


c 


5.05(1 - 

lim G((p) 

^ y 1 - cos 9 


(C22) 


where represents the start of the rolling -up process. 


A simple numerical method for evaluating the denominator of equa- 
tion (C22) is as follows. The loading distribution is given hy the 
series 


G(cp) = 



(C23) 


2 

m+1 


sin sin 
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With equation (C23) ^ the denominator of equation (C22) becomes 


lim G(9) 2 


m m 


° Vl-coscp m+1 


V P V • /lim sin ^ cp \ 

I “"I 


n=i ^^ 1 =^ 


m ^ m 




(C24) 


n=l 


di=l 


For symmetric loading for which and is only odd, 

equation (C24) becomes 


lim G(9) ^ ( 2>/2 

a/i - cos 9 \ m+1 


m 


X >'1=^'' “if Wa±iN * 

li^=odd y \ 2 J 


m -1 

2 


m 


h>/2 

“571 ) ) Gn 


(C25) 


n=i 




=odd 


m 

The summation, ^ sin liiCp^^ in equation (C25) can he evaluated as 


follows : 
let 




i9 
^n 

Z = e = cos cpn i sin 9 


n 


then. 


m 


^ sin n^9n 

|j.i=odd 


imaginary 
part of 


m 




d ) Z 
{T^od-d 
dZ 


-(m+l)(-l) 
2 sin cpj^ 


n 


(C26) 
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With equation (C 26 ), equation (C 25 ) becomes 


m-i 

2 


lim G(cp) 


yiT 


COS cp 


= - G, . - 2 /¥ 


\ 2 / 


n=i 


(-1) Gn 


sin cpj^ 


Equation (C 2 T) inserted into equation (C 22 ) gives 


“ ^o 


5 - 05 ( 1 - 


m~i 

^ r— 'T' (-l)’^^n 

-■/2 <Vm-n-\ - „ 

(^— j sin Vn 


For m = 7 equation (C 28 ) simplifies to the following: 


“ ^o 




, 3/2 


1.4630 Gi - 0.7917 Gg -f 0.6060 Gg - 0.2799 G4 


- ■.-’■'•(4 


0.7315 Ki - 0.3959 K2 -I- 0.3030 K3 - 0.1400 K4 


where the subscripts pertain to span stations- n- 1 , 2 , 3 > ^ 

q = cos — = 0.9239, 0.7071, 0.3827, and 0 , respectively. 

8 


( 027 ) 


(028) 


; or 
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TABLE I.- LATERAL INTERPOLATION FACTORS 



0.098 

0.195 

0.290 

0.556 

0.831 

0.981 

Hi 

H2 

H3 

H4 

H1+H2 

0.0030 

-.0178 

.1086 

.9061 

-.0148 

0.0089 

-.0538 

.3916 

.6533 

-.0449 

0.0102 

-.0645 

.7361 

.3182 

-.0543 

- 0.0449 

.4365 

.8790 

-.2706 

.3916 

0.3378 

.9777 

-.5861 

.2706 

1.6982 

-1.3604 

1.3155 

-.6533 



TABLE II,- VERTICAL INTERPOLATION FACTORS 



±0.1 

±0.2 

±0.3 

±0.8 

X 

±0.1 

±0.2 

±0.3 

±0.8 


For all T] 

n 

“H = 0.707 

Cl 

C2 

0.9611 

.0522 

0.7819 

.2784 

0.5005 

.5986 

-0.1947 

.9978 

1 

2 

3 

4 

1.5928 

-2.7584 

1.4504 

-.2506 

1.7007 

-3.2636 

1.7888 

-.3236 

1.1502 

-2.3921 

1.3447 

-.2566 

-0.4075 

1.0580 

-.6115 

.1384 

n 

T) = 0 

1 

2 

3 

4 

.0973 

-.2957 

1.6725 

- 1.4618 

.0932 

-.2788 

1.9224 

- 1.7412 

.0559 

-.1699 

1.3752 

-1.2816 

-.0029 

.0591 

-.5490 

.5816 

n 

*0 - 0.924 

1 

2 

3 

4 

-11.8776 

6.6138 

-3.5043 

1.5814 

-17.8814 

11.9^^06 

-7.9690 

3.7283 

-15.0373 

10.8359 

-7.7815 

3.6751 

8.4886 

-6.7796 

5.3097 

-2.5254 

n 

Dn, Ti = 0.383 

1 

2 

3 

4 

-.1667 

1.0587 

- 1.8423 

.9066 

-.1392 

1.1846 

-2.1540 

1.0512 

-.0734 

.8303 

-1.5653 

.7570 

.0125 

-.3137 

.6952 

-.8976 
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TABLE III.- FOR SEVERAL VALUES OF m AND M 







|j.( for m 

= 7) 






M = 7 

M =15 

M= 31 

n = 1 

n = 2 

n = 3 

n = 4 




0 


0 

0 

2.6131 

- 1 . 4 i 42 

1.0824 

-0.5000 




— 


— 

1 

4.5889 

-2.1053 

1.5663 

-.7193 







1 

2 

2.8844 

-.2363 

.0751 

-.0253 









3 

.6573 

2.0046 

-1.5402 

.7109 




1 


2 

4 

- 1 . 4 i 42 

3.6955 

- 2 . 4 i 42 

1.0824 







■ 

5 

-2.7625 

4.1506 

-2.0037 

.8553 




— 


3 

6 

-3.1207 

3.1958 

-.3621 

.1005 








7 

-2.5843 

1.2224 

1.8601 

-.8209 




2 


4 

8 

-1.5307 

-1.0000 

3.6955 

- 1 . 4 i 42 









9 

-.4383 

-2.6609 

4.2950 

- 1.2844 







5 

10 

.3114 

-3.2465 

3.3216 

-.3367 








11 

.5651 

-2.7288 

1.1207 

1.1492 




3 


6 

12 

. 4 i 42 

-1.5307 

- 1 . 4 i 42 

2.6131 






— 

13 

.1016 

-.2939 

-3.2260 

3.4397 




— 


7 

14 

-.1258 

.4372 

-3.5579 

3.2212 




— 



15 

-.1444 

.4635 

-2.2904 

1.9416 








m = M = 

= 15 




1 ^ 

n 

= 1 

n = 2 

n = 3 

n = 4 

n = 5 

n = 6 

n = 7 

n = 8 

0 

5.1258 

-2.6131 

1.800 

- 1 . 4 i 42 

1.2027 

-1.0824 

1.0196 

-0.5000 

1 

-2.6131 

6.9258 


■ 4.0273 

3.0027 

-2.4966 

2.2223 

-2.0824 

1.0196 

2 

-3.3258 


l. 4 i 42 

6.3285 

-3.6955 

2.8196 

- 2 . 4 i 42 

2.2223 

-1.0824 

3 

1.1989 


3.9231 


•1.0824 

6.1454 

-3.6131 

2.8196 

-2.4966 

1.2027 

4 


5973 

1.5307 


.4.1062 

-1.0000 

6.1454 

-3.6955 

3.0027 

- 1 . 4 i 42 

5 

i 

3318 


-.7804 

1.6131 

-4.1062 

-1.0824 

6.3285 

- 4.0273 

1.8000 

6 

-.1831 


.4142 

-.7804 

1.5307 

-3.9231 

- 1 . 4 i 42 

6.9258 

-2.6131 

7 

i 

,0824 


-.1831 

03318 

-.5973 

1.1989 

-3.3258 

-2.6131 
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Plan view 
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z or b 


Side view 

Coordinates of flat- sheet vortex system 




Figure 1.- Illustration of coordinates. 
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symbol method 

one lifting line 

three lifting lines 

five lifting lines 






uniform chord loading 


r~ ^j‘£ f sernispans 
distance behind wing trailing edge 



(a) A = 1.0 

Figure 2.- Comparison of the downwash field due to a lifting line with 
that due to a lifting surface composed of several lifting lines. 

(a = 2.0, A = 56°, ft = 0) 
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symbol method 

one lifting line 

three lifting lines 

lifting lines 



cotangent- type chord loading 



uniform chord loading 




T- , semispans 
distance behind wing trailing edge 



(b) A = 0 

Figure 2.- Concluded. 
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(a) T| = 0 

Figure 3*~ Downwash influence coefficient, asn, due "to total flat-sheet vortex system; t/b aft 

of c/4, as=o. 
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(b) n = 0.383 
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Figure 3«- Continued 


MCA TN 33^^6 



4 .6 .8 1.0 

t//3, semispans 


.6 .8 10 
t//8, semispans 


- 1.0 


- 1.2 


- L 4 


- 1.6 


- 1.8 


- 2.0 


MI, VOVM 




.4 .6 .8 1.0 1.2 

T//S, semispans 


14 16 O .2 

(c) = 0.707 

Figure 3»- Continued. 
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Figure 3*- Concluded 
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Figure 4 .- Downwash influence coefficient, S’sn> to total flat-sheet vortex systemj t/p aft 

of c/k, «s = 0.5. 
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(b) Ti = 0.383 
Figure 4.- Continued 
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Figure k.- Concluded. 
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Figure 5.- Downvash influence coefficient^ ^Tn^ 

sheet. 


due to trailing vortex 
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Figure 5-- Continued 
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Figure 5>- Concluded 
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Figure 6.- Downwash parameter for two trailing vortex cores that extend 
downstream from the quarter chord of the wing tip. 
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Figure 6.- Concluded 
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Wing Characteristics 

A‘4 C^=0.5t 

A = 40** a =7.9** 

A = 625 

Longitudinal Position 

2.0 c aft of ^ 
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O experimental downwash 
from contour plots 

—W— downwash computed using 
Weissinger span loading 

~E downwash computed using 

experimental span loading 

•• — experimental wake ^ center 
location 
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span hading 


Figure 7.- Comparison of downwash angles obtained experimentally with those computed hy flat -sheet 

theory. 
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Figure 8.- Comparison of Weissinger span loading and experimental span loading for the configura- 
tion of figure ?• 
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(a) Ti = 0 

Figiire 9.- Comparison of experimentally obtained downwasb angles for a 
ving-fuselage combination with those computed by the flat-sheet 
method plus a fuselage correction. 
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Figure 9*~ Concluded. 
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experiment, references !2, !3 , !4 

computed with sweep correction, eq. 26 

computed without sweep correction, eq. 23, 

assuming the vortex starts at quarter chord 
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Figure 11.- Comparison of computed and experimental spanvise location, 
of the tip vortices behind several wings. 
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(a) A = 2 

Figure 12.- Contoiirs of constant downwash angle, degrees, at 2.0 M.A.C. 
aft of the quarter chord of the wing M.A.C. predicted hy rolling -up 
method and hy flat-sheet method; A = 40°, A = 0.5; A = 2 and 4, and 
Cl = 0.4, and 0.6. 
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Figure 12.- Concluded. 
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